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ABSTRACT 

In this paper we analyze high time resolution single pulse data of PSR 
B0809+74 at 820 MHz. We compare the subpulse phase behavior, undocu- 
mented at 820 MHz, with previously published results. The subpulse period 
changes over time and we measure a subpulse phase jump, when visible, that 
ranges from 95° to 147°. We find a correlation between the subpulse modulation, 
subpulse phase, and orthogonal polarization modes. This variety of complicated 
behavior is not well understood and is not easily explained within the frame- 
w ork of existing models, most of w hich are founded on the drifting spark model 
of iRuderman fc Sutherland! (119751). We quantita tively fit our data with a non- 
radial oscillation model (jClemens fc Rosen! |2008| ) and show that the model can 
accurately reproduce the drifting subpulses, orthogonal polarization modes, sub- 
pulse phase jump, and can explain the correlation between all these features. 



Subject headings: pulsars:individual:PSR B0809-|-74-pulsars:general — pulsars:polarization- 
st ars : neutron — st ars : oscillations 



1. Introduction 



First discovered in 1968 bv !Cole fc Pilkington! (11968!), PSR B0809+74 is a bright, slow 
pulsar with drifting subpulses that has been continuously studied over the past 40 years. The 
literature contains a wide range of behavior including changes in subpulse period, subpulse 
phase behavior, average pulse shape, and orthogonal polarization modes as a function of radio 
frequency. The average pulse profile and polariza t ion an gle histogram for our observations 
at 820 MHz are shown in Figure HJ iHobbs et al.l (|2004[ ) have monitored this pulsar for at 
least 6 years, measuring a spin period (Pi) of 1.292 seconds and a dispersion measure of 
6.116 pc cm~^. A list of basic parameters are in Table H] 
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Fig. 1. — Top: The average total intensity (solid line) and linear polarization (dashed line) 
of PSR B0809+74 from epoch MJD 54922 at 820 MHz, consisting of 232 pulses. Bottom: A 
2D histogram of the polarization angle for the same data. The center of the pulse profile is 
arbitra ry; we chose maximum to b e 60° in pulse longitude to directly compare our data to 
that of Edwards fc Stappers J2003h (see Figure [2]). 
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Property 


Value 


Uncertainty 


PSRJ 


J0814+7429 




RAJ 


08:14:59.50 


0.02 s 


DECJ 


+74:29:05.70 


0.11 " 


Period (-Pi) s 


1.292241446862 


3*10-12 


Dispersion Measure (DM) pm cm~'^ 


6.116 


0.018 


Epoch 


49162.0 




Spin Down (P) s s^^ 


1.68114 X 10-16 


1.4 X 10-20 



Table 1: The basic parameters of PSR B0809+74 faobbs et al.ll2004h 
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Bartel et al. 


1981 


Davies et al. 


1984) 



two adiacent subpulses i n the same pulse and is usually referred to as P2. As discussed in 



Clemens fc RosenI (j2004( ). the value of P2 is not an accurate measurement of the underlying 
fundamental subpulse period, Pume- At higher frequencies, the P2 appears t o decrease to 



39 m s 31 ms, and 29 m s at 406 MHz, 1412 MHz, and 1720 MHz respectively flDavies et al. 



1984J : iBartel et al.l Il98ll). Our measureme nts at 820 MHz fall in the middle range of ob- 



servational frequencies. iBartel et al.l ( 1l98ll ) find that while the subpulse period appears to 
change between 102.5 and 1720 MHz by a factor of 1.8, the time it takes for a subpulse to 
return to the same longitude, P3, remains constant. 

The likely underlying cause for the change in the measured subpulse period is a subpulse 



phase discontinuity (or jump) that appears at h igh frequencies but not at 
The subpu^ 



se phase jump is not seen at 32 8 MHz (lEdwards fc Stappers 



408 MHz (iProszvnski fc Wolszczani 
MHz. lEdwards fc StappersI f l2003h : 



1986), or at 500 MHz f Wolszczan et al. 



2003 



ow frequencies . 



Edwardsll2004[) 



Edwardsl (120041 ) report a phase jump of 



198l|). At 1380 
^ 120° as shown 

in the middle panel of Figure |2J Furthermode, the subpulse phase jump occurs starts at 
approximately 56.5° in pulse longitude (bottom left panel of Figure [2]), corresponding with 
a minimum the subpulse an iplitude envelope (dark line in the top left panel of Figure [2]). As 
Edwards fc StappersI (|2003[ ) discuss, in any given pulse with two subpulses, the subpulses 
generally lie on opposite sides of the subpulse phase jump and appear closer together in pulse 
phase than they normally would in the absence of the phase jump, resulting in a smaller 
value of P2. 



Most drifting subpulse models are based on the drifting spark model ( iRuderman fc Sutherland 



I975I ) in which a vacuum gap forms between the stellar surface and co-rotating magneto- 
sphere due to the charge depletion from the emitted particles. To prevent the vacuum gap 
from growing indefinitely, sparks discharge across the vacuum gap. These sparks are fixed 
relative to each other and form a carousel that rotates around the magnetic pole at a rate 
incommensurate with the spin period of the star. The drifting subpulses are the manifesta- 
tion of these spark discharges. Known as the drifting or (ro tating) spark model, this model 
is the basis for man y current models of drifitng subpulses (jKomesarofI 1970; Backer 1976 
Gil fc Sendvljboool ). 



In I Clemens fc RosenI (l2004l I2OO8I). we propos ed a non- radial oscillation mode l based on 



astero seismological techniques ( jPziemb owskill 1 9 771 ) (see also lRobinson et al.l ( 1l982l ): IClemens et al. 



( 120001 )) as an alternative to the drifting spark model. Pulsations in stars are not uncommon: 
white dwarf stars , ZZ Ceti stars, rapidly oscillating AP stars, and delta Scuti stars a l l show 
pulsation modes (IKleinman et al.lll998[ Ivan Kerkwijk et al.ll2000l ; IWinget et al.lll98ll : iKurtz 
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Fig. 2. — Top left: The average profile and subpulse longitude envelope for PSR B0809+74 
at 328 MHz (dashed lines) and 1380 MHz (solid lines). Top right: The average profile and 
subpulse longitude envelope at 820 MHz taken on MJD 54922. Middle left: The subpulse 
phase, plotted multiple times, spaced 360° apart. The white and dark circles are from data 
collected at 328 MHz and 1380 MHz, respectively. The dotted hne shows the phase slope of 
25°. Middle right: The subpulse phase for our data at 820 MHz; the dotted line shows the 
phase slope of 27.8°. Bottom left: The difference between the subpulse phase and that of 
the phase slope, indicating the magnitude of the phase jump. The phase jump at 1380 MHz 
(dark circles) is 120° and the phase jump at 328 MHz (white circles) is plotted twice with 
an offset of 120°. Bottom right: The difference between the subpulse phase and that of the 
phase slope for our data at 820 MHz, results in a phase jump of 116°. Using an alternative 
method of fitting two lines on either side of the phase jump (middle, right panel) rather than 
subtracting the phase slope results in a phase jump of 145° (see §3.2p . The phase jump at 
820 MHz occurs at the same pulse longitude as the phase jump at 328 MHz and 1380 MHz. 
All the plots in the left panels are reproduced from 



Edwards fc Stapperd ( l2003l ). 
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van Horn 



19821: Bregerlll969t ) . We • were not the first pr oponents of a oscillation model for pulsars; [Gold 



(Il980[ ): IStrohmayerl (|1992[ ) all proposed oscillations as an explanation for 
drifting subpulses. However, these papers did not address the wide range of phenomenology 
seen in pulsars with drifting subpulses. 

In this paper, we analyze high quality single pulse measurements of PSR B0809+74 at 
820 MHz. In §2] we discuss our observations and conduct a detailed analysis of the data in 
^ We then explain our model in ^ and examine the data in the context of our model in 
§3 Finally, in §H1 we discuss how our 820 MHz observations compare the observations at 
other frequencies and explain the single pulse behavior within the context of a non-radial 
oscillation model. 



Observations 



The observations of PSR B0809+74 were taken in the spring of 2009 with the 100-m 
Green Bank Telescope (GBT) using the new pulsar backend GUPPI in filterbank mode. 
The dates and lengths of the 10 observations are listed in Table [21 Full-Stokes spectra were 
acquired in a 200 MHz-wide band centered at 820 MHz radio frequency. The frequency 
resolution was ~98 kHz and the spectrum integration time was 160 /is. The filterbank data 
were then averaged into 1024-bin single-pulse profiles using the ephemeris given in Table [T] 

Flux and pola rization calibration were performed using the PSRCHIVE software package 
(iHotan et al.l 120041 ). Each of the 10 observations was performed at a different hour angle. 
The rotation of the source with respect to the telescope was us ed to solve f o r the receiver 
system's intrinsic polarization cross-coupling matrix, following Ivan StratenI (120041 ). From 
the calibrated profiles, we determined a rotation measure (RM) of —12.2 ± 0.2 ra d m~^, 
consistent with the catalogued value of —11.7 ± 1.3 rad m~^ (IManchesterl Il972l ). The 
calibrated single-pulse profiles were then RM-corrected and integrated over the full band for 
the analysis described in the following sections. 



3. Data Analysis of PSR B0809+74 

3.1. Subpulse Period 

The ten observations we acquired in April 2009 varied in length and morphology, as 
shown in Table [21 Eight of the ten observations displayed nulling, periods with zero emission. 
Only two epochs, MJDs 54922 and 54944b, exhibited a single subpulse period in the FFT 
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of the entire run; the remaining observations had multiple significant periods. The first and 
last observations (MJDs 54922 and 54961) had the brightest flux. When displaying data and 
our fits to the data, we use epoch MJD 54922 as it has the second largest flux and is not 
affected by nulls as is data from MJD 54961. 
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Date 


Epoch 


Length of Observations 


Length of Nulls 


Multiple Periods 


Mean Flux 




MJD 


Pulses 


Pulses 




mJy 


04-01-2009 


54922 


232 


2,2 


no 


310.3 


04-22-2009 


54943 


466 





yes 


91.9 


04-23-2009 


54944a 


397 


3 


yes 


191.9 


04-23-2009 


54944b 


476 


2 


no 


189.9 


04-27-2009 


54948a 


475 


3, 2,2 


yes 


126.5 


04-27-2009 


54948b 


475 


3 


yes 


171.7 


04-27-2009 


54948c 


476 


2, 2, 2, 2 


?? 


114.3 


04-27-2009 


54948d 


476 





yes 


60.2 


04-27-2009 


54948e 


477 


3 


yes 


40.3 


05-10-2009 


54961 


475 


8,4 


yes 


374.6 



Table 2: A summary of our observations at 820 MHz. The first and last observations (MJDs 
54922 and 54961) have the brightest flux. When displaying data and our fits to the data, 
we use epoch MJD 54922 as it has the second largest flux and does not display nulls. 
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We find the fundamental subpulse period, Pume^ to be betwee n 48 and 54 ms, wh icfi 



is consistent witli tfie measured value of P2 by iDavies et al.l (Il984f ) ; iBartel et al.l ( 119811 ) at 
low frequencies without the presence of a subpulse phase jump (see §5.2p . A fast Fourier 
transform (FFT) shows that the subpulse alias peaks at approximately 39 ms. However, the 
FFT calculates the period of the subpulse based on the separation of the peaks, and the 
presence of the subpulse phase jump makes the subpulse peaks artific ially closer than the 



would be in the absence of the phase jump (lEdwards fc Stapperd 120031 ) . 



The average pulse profile at 820 MHz is consistent with those at 328 MHz and 1380 
MHz. As shown in the top panels of Figure [21 the subpulse modulation envelope at 820 MHz 
(dashed line, top right panel) resembles that at 1380 MHz but not at 328 MHz (top left 
panel). This difference is the result of at subpulse phase jump (shown in the bottom panels 
of Figure [2]) at 820 MHz and 1380 MHz; we discuss the subpulse phase behavior in §3.21 

Figure [3] characterizes the subpulse behavior for each observation. The mean flux in 
the top panel of Figure |3] is also described in Table [2l The subpulse period can also be 
characterized by P3, which is the time it takes for a subpulse to return to the same longitude 
after successive spins of the star. P3 is not sensitive to which alias we chose from the 
FFT as the subpulse period. To measure P3, we use the longitude resolved fluctuation 
spectrum (LRFS), a Fourier transform calculated at each spin phase. The value of P3 for 
each observation is shown in the middle panel in Figure [31 The bottom panel displays the 
subpulse phase jump; the method for calculating the subpulse phase jump is discussed in 

The top plots in Figure [H shows the LRFS of two different data sets taken at different 
epochs; the bottom panels of the top plots show a peak at P3 = 14.61 seconds, or ~11.31 Pi. 
The bottom plots of Figure [H shows the corresponding driftband plot, a contour plot of the 
data folded at P3. Figure [31 shows P3, calculated from the LRFS, for all the different epochs. 



3.2. Subpulse Phase and Orthogonal Polarization modes 



Each of the ten data sets show a subpulse phase jump at pulse longitude of 56.5° (offset 
by 3.5° from the maximum in the pulse profile), as shown in the right panels of Figure [2l We 
calculate the subpulse phase using the amplitude and phase of the peak of P3 in the LRFS. 
Figure [21 shows the subpulse modulation envelope, the subpulse phase, and the subpulse 
phase with the nominal phase slope of 28.7° /° removed (bottom r ight p anel). Instead of 
fitting a line to the subpulse phase jump like [Edwards fc Stapperd (120031 ) . we fit two lines: 
one on each side of the phase jump. The solid lines in the right third panel of Figure [21 




54922 



54944 54948 



54961 



Epoch 



Fig. 3. — Top panel: The mean flux as a function of epoch. Middle panel: The period 
P3 at each epoch. Bottom panel: The subpulse phase jump; see §3.21 for our methodology 
on calculating the phase jump. 
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Fig. 4. — Top: The longitude resolved fluctuation spectrum from MJDs 54922 (left) and 
54961 (right), where P3 ~11.31Pi. The left panels show the integrated pulse proflle; the 
maximum on the plot is 1100 mJy. The bottom panels show the power in mJy^/Hz. Bottom: 
The same data set folded at P3. 
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illustrates this fitting. The difference between the two parallel solid lines at pulse longitude 
56.5° is the subpulse phase jump; in Figure[2], the jump is 145.1°. The middle panel of Figure 
E] shows the subpulse phase jump for each epoch and the bottom panel shows the average 
slope (taken from the two lines fit on either side of the jump). The magnitude of the jump, 
present in 7 of the 10 data sets, ranges from 95° to 147°. 

The orthogonal polarization modes have distinctly different behavior on the leading 
and trailing edge of the pulse profile; the behavior changes at 56.5° as well, coincident with 
the subpulse phase jump. As shown in Figure [H the left side of the pulse profile appears 
to be dominated by a single polarization mode; the right side of the pulse profile shows a 
combination of two or thogonal p olariz ation modes. This is consistent with the polarization 



behavior described bv Edwards! fl2004f ) at 1380 MHz. 



Our Non-Radial Oscillation Model 



In previous papers, we developed an oblique pulsator model (lKurtzlll982l ) for pulsars in 
which drifting subpulses are produced by no n-radial oscillations whose pe riods are incom- 
mensurate with the spin period of the pulsar ( IClemens fc Rosenll2004i |2008| ). The non-radial 
modes of our model are aligned to the magnetic axis, so in addition to the drifting time-like 
pulses, our model pr oduces longitude station ary variations caused by nodal lines rotating 
past our line of sight (jClemens fc Rosen! l2004j ). Nodal lines are places of unmodulated emis- 
sion and are described by the zeros in the spherical harmonic in Equation [1] The emission 
on either side of a nodal line is out-of-phase and the subpulse period is defined as rate at 
which a region between two nodal lines emits radio emission. 

The pulsations cause displacements of stellar material which modulate the linearly- 
polarized emission, as might be produced by curvature radiation. These displacements have 
a transverse elec tric field vector that point s towa rd the magnetic pole and follows the rotating 
vector model of !Radhakrishnan &: Cookd ( 11969! ) (see Fig ure 151) ■ This mode of po larization, 
the displacement polarization mode, can be described as (IClemens &: Rosen! !2008! ): 



'^l,m=o{Omag) COs{ujt - IpQ - Ipdelay)) (1) 

where \&i,m=o is a spherical harmonic of high i and m = 0. The variable Omag refers to the 
magnetic co-latitude, because the pulsations in our model are aligned to the magnetic pole. 
The subpulse period, cu, is related to the fundamental subpulse period, Ptime= 27r/w. The 
phase term ifj^ allows for the arbitrary phase of the drifting subpulses. The 'ipdeiay term allows 
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for a time lag between the maximum amplitude of the pulsations and emission maximum; 
for non-adiabatic oscillation s, the thermal maximura can lag the displacements in phase and 
i^deiay allows for this effect ( IClemens fc Rosenll2008[ ). The amplitudes aoopA/ and ai^^p^^ are 
to be fitted to the data, as well as Pume^ i^o, and ipdeiay 

The pulsational displacements and their associated velocities move in the plane of the 
magnetic field. Thus, the induced electric field as a result of the velocities {E = v x 
B) is naturally orthogonal to the Radhakrishnan and Cooke vector (see Figure The 
polarization mode that results from the induced electric field due to the velocities is the 
velocity polarization mode and is expressed as: 



A 



VPM 



it) 



^^vpM^n sm(wt - ^o) 



(2) 



mag 



which incorporates the time derivative and the Qmag derivative of Equation [T], as appropriate 
for horizontal pulsation velocities. This equation is analogous to the Vq in equation three of 
Dziembowskil (119771 ). 



Thermal and field emission from the neutron st ar surface can accel e rate electrons along 
open field lines with the formation of a vacuum gap (jJessner et al.ll200ll ). IStrohmayerl (1l992l ) 
proposed that neutron star oscillations could modulate the radio intensity if greater quantities 
of plasma are injected into the magnetosphere during pulsation n i axima , when local heating 
of the stellar surface is greatest. As discussed in I Clemens fc RosenI (120081 ) . we assume that the 
amplitude of the displacement polarization mode follows surface thermal variations caused by 
non-radial oscillations. This means that for non-adiabatic oscillations, the thermal maximum 
can lag the displacements in phase, parameterized by '^deiay in Equation [H 

The parameters in our model can be divided i nto two groups: geometrical parameters 
and pulsational parameters (IRosen fc Clemensll2008l ). The geometrical parameters, discussed 
in §5.11 are largely independent; only /9, the angle between the magnetic pole and our line- 
of-sight, and £, the degree of the spherical harmonic, are related. There are 7 pulsational 
parameters - amplitudes, periods, and phases and all of these are fit within our model. 



5. Quantitative Fitting of PSR B0809+74 

The fitting process occurs in two steps: we fit the polarization angle swing to determine 
the pulsar geometry and then we fit the data to our non-radial oscillation model to determine 
the pulsational parameters. Since we do not expect the geometry of the star to change over 
time, we added the polarization angle from each of the data sets together determined the 
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geometry from a single fit. Because our model does not incorporate circular polarization, we 
removed the Stokes paramter V from the total intensity such that Imw = -\//^ — V^. 



5.1. Pulsar Geometry 

We fit the polarization angle swing to determine the pulsar geometry, namely a (the 
offset in rotation and magnetic axes), /3 (the angle between the magnetic pole and our line 
of sight), (po (the rotational longitude of the magnetic pole), and Xo (the position angle 
of the linear polarization at (po)- The only pulsational parameter that interacts with these 
parameters is the spherical harmonic degree i because the positions of the nodal lines that 
encircle the magnetic pole are related i, and the path our sightline crosses through these 
nodal lines is related to /3. 

To determine the geometry, we use th e polarization angle histog ram rather than the 



individual polarization angle measurements (lEverett &: Weisbergll200ll ). To do this, we com- 



pute a single histogram of the polarization angle using all ten data sets. We then discard 
all polarization angles less than 90°, so that we fit a single polarization mode. Since the po- 
larization angle is a histogram, each bin has a distinct number of counts. Using the counts 
in each bin (not counting the discarded polarization angles less than 90°), we calculated the 
standard deviation for the counts per bin. We fit only the polarization angle bins in the 
histogram that have counts greater than 2a. The top panel of Figure |6] shows the average 
number of counts per bin for each polarization angle. The dotted line represents our 2a cut 
off; we only fit values above this line. This process prevents our fits from being dominated 
by the noise. 

To fit our data, we convert each count in each bin in the polarization angle histogram 
above the 2a threshold to x-y values. If a bin at a given polarization angle has 20 counts, we 
create 20 x-y pairs at that value. The more counts per bin, the more x-y pairs are created 
at the polarization angle. We do this for all the bins above the threshold and then fit all the 
x-y points. 



In iRosen fc Clemens! (120081 ). we found that the four geometrical parameters (a, /3, Xo, 
0o) cannot be fit independently and it is the ratio of a to P that is significant; for any 
value of a and Xo, corresponding values of /3 and 0o could be found with approximately 
the same goodness of fit. Therefore, to determine the geometry, we fit a and Xo for various 
values of /3 and 0o- Figure [7] shows a map of the standard error for all trial values of /3 and 
(po- We determined the best geometry from the combination of these two parameters that 
has the smallest standard error, denoted by the circle in Figure [3 
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We find that a /3 of —3.1° and a (j)o of 76.75° liad tlie smallest standard error. For 
these values of (3 and (po , the best fit results in an a of 16.81° and a Xo of —199.29°. 
Interestingly, the fit places the closest crossing of the magnetic pole, (po, on the trailing 
edge of the pulse profile, offset from the maximum intensity. This is similar B0656+14: 
the average polar ization angle shows a shal l ow slope slightly mor e curv ed at the tail end of 
the pulse profile (lEverett fc Weisbergll200ll ). Everett fc Weisbergj (120011 ) fit the polarization 
angle and find 0o to be offset from the maximum in the total intensity by 14.9° in pulse 
longitude; in our case, 0o is offset from the maximum in the total intensity by 19.17°. 



These v alues dif f er frorn that of iRankinI (Il993bl ). who found a and /3 to be 9° and 4.5°, 
respectively. iRankinI (jl993a|) determines the geometry based on the half-power widt h of the 
pulse p rofiles and the beam radius at 1 GHz. The difference between the values of [Rankin 
(Il993bl ) and our fits is most likely due the different approach to det ermining the geom etry 
and to the change in geometry with respect to observing frequency (ISmits et al.ll2006l ). We 



find that the value of 0o, while at the trailing edge of the p ulse profile, is consistent wit h 



the standard polar cap size. Using similar analysis to that of IWeltevrede fc Wrightl (120091 ) . 
we find that the emission height is about O.OSRlc, where Rlc is the light cylinder radius, 
and that the size of the magnetic cap is between the magnetic axis and the last open field 
line. Even if we assume that the steepest polarization angle swing coincides with the center 
of the pulse profile and half of the emission region is missing, the polar cap size does not 
extend past the last open field line. 

In the bottom panel of Figure O we show the portion of the histogram data that we 
use to fit polarization an gle as well as our fit. We also show the polarization angle using the 
values in IRankinI (Il993bl ) and assuming 0o is at the maximum of the pulse profile (~ 18° in 
Fi gure |7| o r 56.5° i n pulse longitude) and Xo = —134°. However, for th e fit using the values 



m 



RankinI (Il993bl ). we fiipped the sign on f3 for internal consistency (lEverett fc Weisberg 



20011). 



5.2. Pulsational Parameters 

To verify that a longer subpulse period, consistent with the literature at low frequencies, 
matches the data better than a period of 39 ms, we fit our model to the data which we discuss 
in detail in §5.21 For each of the 10 data sets, we fit the full range of observed subpulse aliases 
from 31.5 to 51.5 ms. For 9 out of 10 scans, a smaller root mean square (RMS) residual is 
obtained for a subpulse period in the range 48 to 54 ms than the 39 ms where the power 
in the FFT peaks. We fix the subpulse period in our fits rather than letting it be a free 
parameter because several data sets show that the peak due to the subpulse period has a 
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secondary peak of lesser amplitude at a spacing several time larger than the resolution in the 
FFT. We speculate that the presence of nulls in some of t he data sets causes the sub pulse 



period to wander throughout the observations (see Table [2]) ■ Ivan Leeuwen et al.l (120031 ) show 



similar behavior where the subpulse drift rate in PSR B0809+74 changes after nulling. 

After establishing the pulsar geometry, we fit the pulsational parameters, namely: the 
spherical degree, i, the amplitudes of the displacement and velocity polarization modes 
i^iDPAn ^OopAn oo^p^j^), the arbitrary phase of the drifting subpulses (ipo), and the phase 
that allows for a time lag between the maximum amplitude of the pulsations and emission 
maximum (tpdeiay)- Because the subpulse period has has a small secondary peak in some of 
the data sets (see §3.ip . we fix the subpulse period based on the subpulse alias in the FFT 
(see below). The subpulse period is similar but unique for each data set. 

As with our fits to the polarization angle, we fit each data set to our model for various 
values of i. For each data set, the value of i that returns the smallest is the best fit 
to the data. We treat each data set separately in determining the goodness of fit because 
of the variation of stochastic pulse amplitudes. The data sets that have small variations in 
pulse amplitude will have a smaller standard errors compared to data sets that have large 
variations. We estimate the noise, CTj in our data using the radiometer equation: 



sys 



(3) 



where Tgys is the system equivalent flux density (SEFD) in mJy (~ 14 — 18 Jy depending on 
the epoch), B is the bandwidth (200 MHz), and r is the integration time (160 /is). 

The top panel Figure [8] shows total for all the data; the bottom panel shows the 
individual Ax^ for each epoch as a function of spherical degree, i. Each dot represents one 
trial, integer value of £. For all data sets, i is either 18 or 19; the difference is within the 
error of our fit. The best value of £ based on the total is 19. 

We use a method similar to that for finding the best value of i to find the best subpulse 
period. Since our model fits the data in real space (J, Q, and U) rather than Fourier space, 
and a subpulse phase jump is a natural part of the model due to the presence of a nodal 
line, fitting the data to our model is more accurate method for determining the subpulse 
period than using the FFT. For each data set, we choose a subpulse alias from the FFT and 
fit our model to the data, calculating x^. We repeat this for every subpulse alias from 31.5 
to 58.5 ms for each data set, using the value of the subpulse period with the smallest as 
the best fit. The top panel of Figure [H] shows the total for values of the subpulse period 
using different subpulse aliases in the FFT from 31.5 to 51.5 ms; each dot is one alias from 
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the FFT. The bottom panel shows the individual Ax^ for each epoch. The subpulse period 
for all data sets except for one is between 48 and 54 ms; the best subpulse period based on 
the total is 51.5 ms. 

Once we determine the best value of i and the subpulse period for each data set, we 
use Equations [1] and [2] to fit our single pulse data to the model. From Equations [1] and 
[21 we can create the primed Stokes param eters, the Stokes parameters in the non-rotating 
frame of the star jRosen fc Clemens 2008 ). However, instead of defining an emission window 
as we did in iRosen fc Clemens! (120081 ) . we multiply t he Stokes parameter / b y the average 
pulse shape. We then used equations 7, 8, and 9 from lRosen &: Clemens! (120081 ) to transform 
the primed Stokes parameters into the observer's frame of reference (unprimed space), and 
similarly multiplied L, which is invariant under this transformation, by the average linear 
polarization. For each epoch, we use the best value of i and the subpulse period, as 
determined in the process outlined above, and fit the remaining parameters (ciiopa/j '^Oophn 
^OvpM^ V'o, 8i.iadip delay)- Figure [TOl shows the reduced ix^u) fo^ ^^^h epoch. The value of 
x'^i^ strongly correlates with the mean fiux (see the top panel of Figure [3]). Table [3] lists the 
best fit for the geometrical and pulsational parameters for MJD 54922. 



Fig. 5. — The polarization geometry near the surface of a neutron star. The magnetic field 
Bstar extends outward, normal to the stellar surface. The electric field has two components: 
points in a longitudinal direction and is oriented in a latitudinal direction with respect 
to the magnetic pole. The dominant velocity vector v points in the ± 9mag direction, toward 
and away from the magnetic pole. 
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Fig. 6. — Top panel: The average number of counts per bin for each polarization angle in 
the histogram. The dotted line represents the 2a threshold; polarization angles below this 
threshold are not fit. Bottom panel: The fitted portion of the polarization angle histogram 
(above the threshold) using the combined data from all ten data sets. The solid black line 
is the result of our fitting process; the dashed line is using published values (see text). 
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Fig. 7. — The grayscale contours show a for each of the values of /3 and 0^. The solid 
lines trace the standard error. We used increments of 0.1 and 0.25 degrees for /3 and 
0o respectively. The black circle marks the values of /3 abd 0o which result in the smallest 
standard error. 
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Fig. 8. — Top: The total plot of the for various values of £, in increments of one. Bottom: 
A plot of the Ax^ for each epoch. 
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Fig. 9. — The total plot of the for various values of the subpulse period. Each dot 
represents one alias of the subpulse period taken from the FFT. We fit subpulse aliases from 
31.5 to 58.5 ms. Bottom: A plot of the Ax^ for each epoch. 
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Fig. 10. — We fit the data to our model, where the amphtudes and phases are the free 
parameters. We set the values of ^ and the subpulse period based on the methodology 
described in the text. The plot shows reduced hdv) resulting from out fit for each epoch. 
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Parameter 


Value 


a 


16.81° 




-3.1° 


(f>o 


76.75° 


Xo 


-199.29° 




1.201 




0.733 


^OVPM 


0.055 


V'O 


-85.66° 




113.857° 


P2 


51.503 ms 


I 


19 


X 1/ 


7.24 



Table 3: The geometrical and pulsational parameters resulting from our fit of the data from 
54922 to the model. 
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Using the fitted parameters, we create synthetic hghtcurves and compare them to the 
data. Figure [11] shows data from MJD 54922 and a simulation using our fitted parameters. 
The data shows a phase jump of 145.1° while the data produced by our model has a phase 
jump of 187.7°. The fact that our model does not reproduce the phase jump exactly is not 
surprising; the phase jump is fit indirectly in our model. Our model fits /, Q, and U directly 
and the phase is computed secondarily. The phase jump occurs at 56.5°, coincident with the 
nodal line, as a result of our best fit of £. Figure [12] shows the longitude resolved fluctuation 
spectrum and driftbands of the data from MJD 54922 and simulation. The harmonic in the 
simulation of 0.14 Hz is more pronounced than that of the data. This is mainly due to the 
fact that the data neither contains pure sinusoids nor a perfect window function, while our 
simulation does. The two phase terms in our fitting, ipo and ip delay, are most prominent in 
the driftband plot. In the data, the maximum intensity occurs between a pulse phase of 60° 
and 65°; altering ipo changes the location of the maximum intensity in both pulse phase and 
P3. In the simulation, the driftband has a feature at a pulse phase of 65° and period of 14 
seconds. This small feature, less obvious in the data, is due to the velocity polarization mode. 
The value of ipdeiay dictates the spacing of between the two polarization modes. Figure [T3] 
shows single pulses of the same data and simulations as in Figure [TT] In the simulations, the 
subpulses on the leading edge of the profile are due to the velocity polarization mode and 
are slightly weaker than those due to the displacement polarization mode. 



6. Discussion 



6.1. Observations In the Context of a Non-radial Oscillation Model 



A non-radial oscillation model explains the wide range of behaviors seen in slow pulsars 
and is substantially simpler and more cohesive model than those based on drifting sparks. 
We expect that the subpulse period should not change with radio frequency. We also expect 
that a nodal line must correspond with a region of zero modulated emission and a subpulse 
phase jump. 



The radio frequency dependence of is a matter of debate. jPavies et al.l ( 119841 ) found 
that for PSR B0809+74, th e subp ulse spacing (P2) appears to change with frequency, but 
P3 does not. Ilzvekova et al.l (119931 ) found similar result s for four additional p ulsars. However, 
the average pulse shape also changes with frequency (Ilzvekova et al.lll993l ). resulting from 
the divergence of the magnetic field and the change in radial distance from the magnetic 
pole. High frequencies originate close to the stellar surface, so the average pulse profile is 
narrow, but emission at these frequencies originates further from the rn agnetic pole radially 



than low frequencies (iKomesarofi 



I970I : ICordeslll978l : ISmits et all 120061 ) 
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Fig. 11. — Top panels: The average pulse profile (solid line), the average linear polarization 
(dashed-dotted line), and the subpulse modulation envelope (dotted line) of the data from 
MJD 54922 (left) and model (right). Middle panels: The polarization angle histogram for 
both the data and model. Bottom panel: The subpulse phase, calculated from the LRFS 
(see §3.2p . The solid line is a linear fit to both components of the subpulse phase. 
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Fig. 12.— Left column: the LRFS (top) and driftband of data from MJD 54922. The left 
panel shows the integrated pulse profile; the maximum on the plot is 1100 mJy. The bottom 
panel shows the power in mJy^/Hz. Right column: the LRFS and driftband of simulated 
data using the best parameters as determined by our fit of the data to our model. The panels 
in the LRFS have the same scale and units as the data in the top left plot. 
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Fig. 13. — Left panels: A sequence of 100 pulses from data from MJD 54922 (lower panel) 
and the average profile (upper panel). Right panels: A sequence of 100 pulses created from 
the fitted parameters from our model against the data from this epoch. 
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This scaling law, called radius-to-frequency mapping, can be mistakenly applied to 
in two ways. As the average puls e profile broadens at lower frequencies, the spacing 



Wolszczan et al. 



between the subpulse increases as well (Ruderman &: Sutherland! Il975t llzvekova et al 



198ll : lBartel et alJl98ll : lGil fc Krawczvk 



1993 



19961 ). However, as Edwards &: Stappers 



fc002h describe, the change in P2 with frequency is a result of the change in window function: 
as the average pulse profile broadens (or narrows) the window function that modulates the 
subpulse period changes as well, affecting the apparent spacing between subpulses. Secondly, 
as illustrated in PSR B0809+74, the average pulse profile and subpulse phase behavior can 
change quite noticeably between different frequency. The appearance of a subpulse phase 
jump will change the measured value of P2 (see ^3.11). Both the subpuls e period and phase 
envel ope should be independent of frequency ( lEdwards fc Stappersll2002l : IClemens fc Rosen 
20o3). 



Because the oscillation mode is a fundamental property of the star, we expect the 
subpulse period to be independent of observing frequency. However, the subpulse period 
can be distorted by both the presence of a nodal line and the width of the average profile 
( IClemens fc RosenI |2004J ) , and therefore the apparent spacing of P2 is not an accurate 
measurement of a stable, underlying clock. In the case of PSR B08094-74, si nce the change 
in frequency effectively changes our sightline traverse, (3 ( ISmits et al.ll2006l ). the subpulse 
phase jump that appears at higher frequencies can be attributed to a nodal line moving 
into our line of sight due to the change in effective geometry. The subpulse period should 
remain the same, but the apparent spacing changes. We see this in the period of PSR 
B0809+74: the FFT incorrectly determines subpulse spacing to be around 39 ms, but fitting 
the single pulses more accurately measures of the subpulse period to be 48 to 54 ms. While 
our me asurements of the subpulse period appear to be consistent with p revious pubhshed 
results ( jBartel et al.l Il98ll : iDavies et al.l Il984j : Edwards &: StappersI l2003l ) , the true test to 
determine if the subpulse period is independent of frequency is to conduct simultaneous, 
multifrequency observations. 



The subpulse phase jum p appears to change with observing frequency. Edwardsl (120041 ) ; 
Edwards fc StappersI (120031 ) reports a phase jump of approxima tely 45° and 120° in the tota l 
intensity at 328 MHz and 1380 MHz. Using the methodology of Edwards fc StappersI (120031 ). 
we find the subpulse phase jump to be 116° at 820 MHz; using a slightly different technique 
described in §3.2[ we find the phase jump to be 145°. Regardless, the appearance of the 
phase jump at 820 MHz and above is consistent with a nodal line moving into our line-of- 
sight. This hypothesis is borne out in the subpulse modulation envelope shown in the left 
upper panel of Figure [H The dotted lines in the top left panel show the average profile and 
subpulse modulation envelope at 328 MHz; the solid lines are the same for 1380 MHz. The 
subpulse modulation envelope at 328 MHz is a single peak, while the modulation envelop at 
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1380 MHz shows a minimum which would be consistent with the presence of a nodal line. 



In our model, the choice of I =19 places a nodal line at 56.5° and produces the phase 
jump at the same location. The combination of a nodal line and (po offset from the center 
of the pulse profile is responsible for the presence of two orthogonal polarization modes on 
the trailing edge of the pulse profile. While the amplitude of the displacement and velocity 
polarization modes are constant throughout the simulation, the asymmetry due to I and 

00 cause the velocity polarization mode to be suppressed on the leading edge. 

This value of ^ is significant ly sr aaller than the value we c alculated for PSR B0943+10 
=75) (iRosen fc Clemens! |2008| ). In I Clemens &: RosenI (120041 ) we estimated the magnifica- 
tion factor to map the £ we observe to the surface of the star. We found that at 1 GHz, 
the magnification factor was ~ 7, and thus in our qualita t ive rn odel an £ of 85 translated to 
an I of 600 on the stell ar surface. In iRosen &: ClemensI (120081 ) . we found I to be 75 using 
430 MHz archival data (IPeshpande fc RankinlboOlh . Assuming that observations at lower 
frequencies are at higher altitudes from the surface of the star, and thus the magnetic field 
is more divergent, the magnification factor is then a lower limit and the C. of 75 translates 
to a lower limit of d. =525. In these data, with an i of 18 or 19, the lower limit becomes an 

1 of 130 on the stellar surface. We note here that the value of £ =75 for PSR B0943+10 
was not rigorously determined using the methodology outlined in this paper, and smaller 
values for I for that star are possible. 



7. Conclusions 

This paper shows t he second quantitative f it of our non-radial oscillation model to single 



pulse data of a pulsar ( IRosen fc ClemensI |2008| ) . In this paper, we 



show the subpulse period and subpulse phase jump, previously unpublished at 820 
MHz. 

show the subpulse phase jump and subpulse period changes with epoch. 

are able to quantitatively determine the best value of £ and the subpulse period. Our 
method for determining the subpulse period is more accurate than using the FFT as 
it determines the subpulse period in real space rather than frequency space, which can 
be affected by the subpulse phase jump. Our fits using J, Q, and U account for the 
subpulse phase jump since it is a natural part of the model. 

are able to quantitatively fit single pulse data to our model and determine a goodness 
of fit using and statistics. 
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• can create simulations based on our fitted parameters which accurately reproduce the 
subpulse period, subpulse phase jump, and orthogonal polarization modes. 

The morphology of PSR B0809+74 is explained easily and naturally within a pulsational 
model. Our non-radial oscillation model is based on established ast eroseismological pri nciples 



that have explained white dwarf variations for the past 40 years (jPziembowskil 119771 ) . This 
model is a viable alternative to the drifting spark model and can provide physical insight in 
the emission mechanism and physical structure of neutron stars. 

The National Radio Astronomy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated Universities, Inc. We would like 
to thank Geoffrey Wright for his thoughtful comments. 
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